Long-time Low-latency Quantum Memory by Dynamical Decoupling 
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Practical quantum memory applications demand high-fidelity storage of arbitrary states for long 
times, while ensuring a small access latency and operating under technological control limitations. 
Building on exact results available for the paradigmatic case of a single qubit exposed to phase 
noise, we identify periodic repetition of a high-order dynamical decoupling sequence as a systematic 
strategy to meet these demands. We demonstrate how numerical search can be invoked to efficiently 
find the sequence best suited for repetition as well as for validating our analytical insights. 
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The preservation of arbitrary quantum states - that is, 
quantum memory - in realistic, noisy physical systems is 
vital for quantum-enabled applications including secure 
communications, quantum computation and simulation, 
and quantum sensing and metrology. While numerous 
techniques relying on both open- and closed-loop control 
schemes have been devised in order to address this chal- 
lenge, quantum error suppression strategies based on dy- 
namical decoupling (DD) [1-4] are emerging as a method 
of choice for physical-layer dccohcrcncc control in non- 
Markovian open quantum systems. In particular, over 
the last few years, different multipulsc DD protocols have 
been investigated in experimental qubit platforms as di- 
verse as trapped ions [5], atomic ensembles [6], semicon- 
ductor quantum dots [7, 8], impurity spins in silicon [9], 
and nitrogen- vacancy centers in diamond [10]. 

The above investigations have consistently pointed to 
DD as a resource-efficient approach to substantially re- 
ducing physical error rates. However, they have largely 
focused on two control regimes [11]: the "coherence- 
time" regime, where the goal is to extend the characteris- 
tic ("1/e" or T2) decay time for coherence as long as pos- 
sible, and the "high-fidelity" regime, where the goal is to 
suppress errors as low as possible for storage times short 
compared to Ti (for instance, during a single quantum 
gate). Still, these two regimes do not capture the typ- 
ical operating conditions of any true quantum memory, 
namely, error probabilities deep below the fault-tolerance 
threshold maintained for arbitrarily long storage times. 
This would be required, for instance, in a quantum re- 
peater, or in a quantum computer where some quantum 
information must be maintained while large blocks of an 
algorithm arc carried out on other qubits. 

In this Letter, we show how DD can be employed to 
realize a useful quantum memory capable of preserving 
quantum information for long times with high fidelity, 
while meeting the practical requirement of a small access 
latency. We consider the simplest yet experimentally rel- 
evant setting of a single qubit undergoing phase deco- 
herence and controlled through sequences of (ideally) in- 
stantaneous 7r pulses. We identify the periodic repetition 



of a high-order DD sequence as an effective strategy for 
memory applications, and analytically characterize the 
resulting error rate and long-time coherence. This is ac- 
complished by working within a filter-design framework 
which directly generalizes the transfer-function approach 
widely used across the engineering community [12] and 
provides a transparent picture of the controlled dynamics 
in the frequency domain, especially attractive for exper- 
imental implementations [5, 11, 13]. Through this ap- 
proach, we identify conditions under which a "coherence 
plateau" can be engineered, and qubit fidelity guaranteed 
to a desired level at long storage times. Finally, we use 
numerical search within the practically important family 
of digital DD sequences to both find optimal sequences 
in the long-time limit, and to verify that a periodic struc- 
ture naturally emerges in the best performing sequences 
at sufficiently long storage times. 

DD for a single qubit undergoing phase decohcrence 
due to the interaction with a quantum bath or, in the 
semiclassical regime, a stochastic field along the (z) quan- 
tization axis, has been extensively discussed in the liter- 
ature [1, 4, 14-18]. Thus, we only summarize key con- 
cepts and notation. Each control operation effects an 
instantaneous tt rotation around the x axis, and the DD 
sequence may be specified in terms of the pulse timings 
{ij}™ =1 , where we also define to = 0, and t n +i = T p 
as the sequence duration. We take all the interpulse in- 
tervals (tj+i — tj) to be lower-bounded by a minimum 
interval t [171. DD implements the control propagator 



U c {t) = aS 



[y P (t) + l]/2 



, where the sequence propagator y p (t) 



is a piecewise-constant bi-valued function that switches 
between ±1 whenever a pulse is applied [19]. 

The effect of DD on qubit dephasing may be evalu- 
ated exactly via calculation of a spectral overlap of the 
control modulation and the noise power spectral density, 
S(lu) = S(—lu), which is determined by the Fourier trans- 
form of the two-time noise correlation function. This is 
true so long as the environment can be approximated 
as a bosonic bath in thermal equilibrium or the fluctua- 
tions are a stationary Gaussian process. In most physi- 
cal scenarios of interest, S(u>) has a power-law behavior 
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at low frequencies, and decays to zero beyond an upper 
cutoff uj c > 0, that is, S(uj) oc \uj\ s J(\uj\,uj c ), and the 
"rolloff function" / specifics the high-frequency behav- 
ior, / = 0(|cd| — oj c ) corresponding to a "hard" cutoff. 
Let y p (uj) denote the Fourier transform of y p (t), which 
can be written as y p (uj) = co^ 1 2?=o(~ I)- 7 '[exp(itju) — 



The filter function of the sequence p 



is then given by F p (uj) 



2 u~, 



y p (uj)\ , and the controlled 



qubit coherence decays over time as e Xp , where the de- 
coupling error x P at time t = T p is 



Xp 



2ituj 2 



F p (lu) dw. 



(1) 



The case n = recovers free evolution over [0,T p ]. For 
any sequence p, the order of error suppression a p is de- 
termined by the scaling of the filter function near u> = 0, 
that is, F p (oj) oc (cjt) 2 ( q p +1 ). A high multiplicity of the 
zero at u> = ultimately leads to a perturbatively small 
value of Xp as l° n g as w c t <C 1. 

Achieving low-error over a desired storage time T s 
is straightforward in principle if arbitrarily fast control 
(t — > 0) is allowed, since it suffices to use a high-order 
DD sequence such as Concatenated DD (CDD, [3]) or 
Uhrig DD (UDD, [4]) over T s = T p . For instance, the 
UDD„ sequence achieves a p = n with n pulses, applied 
1 [7rj'/(2n + 2)]. Increasing a p implies in- 



at tj 



T p sin 



creasing n, at the expenses of shrinking the minimum 
control time scale as r = t\ = 0{T p /n 2 ) [17]. 

Moving from this idealized scenario to a practi- 
cally useful, long-time quantum memory requires several 
linked issues to be considered: (i) Perturbative DD se- 
quences are not generally viable for high-fidelity long- 
time storage. Once the constraint r > is enforced, 
increasing a p necessitates extension of T p . As T p grows, 
however, the underlying perturbative corrections even- 
tually catch up, and a maximal duration emerges be- 
yond which further increasing a p will increase the to- 
tal error [17, 20]. (ii) Numerical DD approaches such as 
randomized [21] or optimized Band width- Adapted DD 
(BADD [17]), which can in principle provide low error 
at long times, become impractical as the configuration 
space of all possible DD sequences over which to search 
grows exponentially with T s . Simplified sequence gener- 
ation approaches are therefore required, (iii) A stored 
quantum state must be accessible not just at a desig- 
nated retrieval time but throughout a time span during 
which it might serve as an input to a quantum protocol. 
Since, however, DD exploits interference pathways be- 
tween control-modulated trajectories, mid-sequence in- 
terruption (t < T p ) may result in significantly sub- 
optimal performance (Fig. 1). Reducing access-latency 
well below the maximum storage time T s , while main- 
taining low error rates, must therefore be considered in 
designing DD approaches to quantum memory, (iv) In 
real settings, the practical challenges posed by sequenc- 
ing complexity should be minimized while maintaining 




FIG. 1: (Color online) Decoupling error and coherence (in- 
set) during a UDD5 sequence with minimum interpulse time 
r = 1 /j,s. Pulse times are marked with filled circles while 
the open circle indicates the readout time T p . Minimal er- 
ror (maximal coherence) is reached only at the conclusion 
of the sequence, with the coherence spike near 2 fis resulting 
from a spin-echo effect. Throughout this work, we assume 
a phenomenological noise model appropriate for a singlet- 
triplet spin qubit in GaAs, S(ui) = auj c (ui /uj c )~ 2 e~ u ' 
with uj £ [wniB,w m ax]. We set a — 0.207, lj c /2it = 10kHz, 
w m in/27r = 0.01Hz, and w max /27T = 10 8 Hz to maximize agree- 
ment with the measured T2 (~ 35 ns) [7, 22] . We chose r well 
above technological constraints (~ ns) in order to reduce n. 



performance. Since the minimum time scale r typically 
is set by the clock speed of the electronic hardware, this 
mandates digital modulation, whereby all interpulse in- 
tervals are constrained to be integer multiples of r [18]. 

Addressing all such issues simultaneously motivates a 
systematic "modular" approach to generate low-error DD 
sequences for long-time storage out of shorter blocks. 
We identify periodic repetition of a base, high-order DD 
cycle as a natural strategy towards meeting our goal. 
In fact, this general approach is widely used implic- 
itly, e.g., the well-known Carr-Purcell (CP) sequence 
{t, 2t, 2t, • • ■ , 2t, t} is generated by repeating the se- 
quence {t, 2t, r}. In order to develop an analytical un- 
derstanding of the effects of repeating an arbitrary DD 
sequence, we observe that filter functions have simple 
transformation properties under sequence combination. 
Imagine that two sequences, p\ and P2, are joined to 
form a longer one, denoted p\ + P2, with a propagator 
y Pl+P2 (t) that is made by continuously attaching each 
sequence propagator. In the Fourier space we have 



yp2^)- 



(2) 



Consider now repeating a sequence p, of duration T p , 
m times, and let [p] m denote the resulting sequence of 
duration T s = mT p . By iterating Eq. (2), the following 
exact expression is found for the DD error: 



X[pV 



S(u}) sin 2 (mujTp/ '2) 



2ttw 2 sin"(u/r p /2) 



F p {uj) duj. 



(3) 



Eq. (3) describes dephasing dynamics under arbitrary 
periodic multipulse control, generalizing the special case 
of Periodic DD (PDD), p = {t, t} [16]: the single-cycle 
filter function F p (lj) is multiplied by a factor which is 
rapidly oscillating for large m and develops peaks scaling 
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FIG. 2: (Color online) DD error for repetition of PDD(o), 
CP(o), and CDD3(n) sequences (a p = 1,2,3, respectively). 
S(cj) oc u) 8 , with s = —2 (solid) and s = (dashed). The 
minimum interpulse interval r is indicated in each panel. 
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FIG. 3: (Color online) DD error for repetition of UDD se- 
quences. The sequences are labeled by both their pulse num- 
ber and duration. The solid series depicts the single cycle 
errors whilst the dashed lines represent sequences repeated 
2, 4, 8, and 16 times (with increasingly dark shading). The 
minimum-error UDD sequences for multiple repeats vs. single 
cycle are distinct, i.e., UDD 5 vs. UDD 7 (UDD 3 vs. UDD 4 ) 
in the left (right) panel. 



with 0(m 2 ) at multiples of the "resonance frequency" 
w Tes = 2ir /T p , introduced by the periodic modulation. 

After many repeats, the DD error is determined by 
the interplay between the order of error suppression of 
the base sequence (F p (lj) oc w 2 ( q p +1 )) ; the noise power 
behavior at low frequencies {S(lj) tx w s ), and the size of 
noise contributions at the resonance frequencies. The sit- 
uation is simplest in the case of a hard upper frequency 
cutoff at uj c . Using the Riemann-Lebesgue lemma, the 
oscillating factor may be removed, resulting in the fol- 
lowing asymptotic expression: 

f bm XW -=XW-=y o w sin 2 W2) ^ (4) 

provided that X[p]°° is finite. To visualize, the DD error 
initially increases as (to 2 x p ), until coherence saturates 
to a non-zero residual plateau at e _x t p )°° , and no further 
decoherence occurs [23]. 

The emeregence of coherence plateau after repeating 
DD sequences can guarantee high fidelity throughout 
long storage times with a latency tt — T p <C T s . Since, 
by comparing Eqs. (3)-(4), it follows that X[p] m < ^X{p\°^ 
for all to, it suffices to pick a base sequence with a suf- 
ficiently small residual error X[p]°°- Mathematically, the 
existence of such a coherence guarantee requires that the 
following conditions be simultaneously obeyed: 

s + 2a p > 1, T p u! c < 2tt, (5) 

which corresponds, respectively, to removing the possible 
singularity of the integrand at or uj Tes in Eq. (4). Thus, 
judicious selection of a base sequence, fixing a p and T p , 
can guarantee saturation of coherence in principle. 

In a realistic scenario where the noise spectrum does 
not have a hard cutoff (e.g., where a residual noise-floor 
persists beyond U) c ), it is impossible to fully avoid the sin- 
gular behavior introduced by the periodic modulation at 
high frequencies as to — > oo. Contributions from the res- 
onating region uj m co rcs arc amplified with to, causing the 



error to increase unboundedly with time and coherence 
to ultimately decay to zero. Nonetheless, a very large 
number to of repetitions may still be applied before the 
resonating contributions become important. This makes 
it possible to engineer a coherence plateau over an in- 
termediate time regime which can be very long from a 
practical standpoint, depending on how fast S(u>) decays 
beyond oj c . Fig. 2 illustrates the emergence of this co- 
herence plateau based on the order of error suppression, 
dp, and noise power law s, as well as (lower panel) the 
exit from the plateau due to the relatively soft Gaussian 
rolloff as a larger r (and subsequently T p ) is used and the 
condition set by Eq. (5) is no longer met. 

Finding the optimal base-sequence for repetition can 
still be challenging and non-intuitive. For example, the 
sequences best suited for long-time storage need not co- 
incide with those optimizing short-time performance. In 
Fig. 3 we compare calculated error rates for single cy- 
cle and multiple repetitions of UDD sequences. The 
local minima indicate that the sequence with minimal 
error at a single repeat is different from the sequence 
with minimal error at many repetitions. In general, nu- 
merical search is necessary to find the optimal base se- 
quence for repetition [17], with a sufficiently small la- 
tency, ti = T p <C T s , to meet system constraints. A 
main advantage compared with a direct optimization over 
all possible sequences running up to T s (BADD) is that 
searching for a low-error sequence with a duration con- 
strained to 2tt/u> c , that will be repeated up to T s , is far 
more tractable computationally. 

Remarkably, a periodic structure for low-error long- 
time DD still naturally emerges through a direct BADD 
search up to time T s , when additional sequencing con- 
straints are enforced. We turn to the Walsh DD (WDD) 
family of sequences to employ a digital modulation that 
is efficient in the complexity of sequencing [18]. Starting 
with a free evolution of duration r, all possible WDD 
sequences can be recursively built out of simpler WDD 
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FIG. 4: (Color online) Top: Minimal-error DD sequences from 
numerical search over WDD sequences, for r = 0.1, 1, 10/xs. In 
each series, the minimal-error sequences systematically access 
higher orders of error cancellation (via concatenation) over 
increasing running times, until an optimal concatenated se- 
quence is found which is then repeated in the longer minimal- 
error sequences. Note that spin-echo, CPMG, CDD and its 
relatives appear naturally in the WDD basis [18]. The gradual 
increase in error for the series with r = 10/is (loss of coher- 
ence plateau) is due to the softness of the high-frequency cut- 
off. Bottom: Control propagators corresponding to the solid 
markers in the middle data series (r = lfis), showing the 
emergence of a periodic structure for sufficiently long storage 
time. Control propagators scaled to same length for ease of 
comparison. Dashed box highlights base sequence CDD4 that 
is repeated for long times. 



point to compensate for finite- width errors, whereas the 
sensitivity of DD to systematic pulse imperfections can 
be largely eliminated by incorporating robust composite 
pulses or symmetric self-correcting cycles, as successfully 
demonstrated in recent DD nuclear magnetic resonance 
experiments with long pulse trains [24]. 

In summary, we have addressed a timely and impor- 
tant problem in quantum information - determining a 
means to effectively produce a practically useful quan- 
tum memory via DD. We have identified the key require- 
ments towards this end, and developed a strategy for 
sequence construction based on repetition of high-order 
base sequences. Our results allow analytical bounding of 
the long-time error rates and identify conditions in which 
a maximum error rate can be guaranteed for long times. 
Finally, we have validated these insights and analytic cal- 
culations using an efficient search over WDD sequences 
under realistic noise conditions. Based on these results, 
we expect that repetition of suitably designed DD cy- 
cles will provide a practical avenue to high-fidelity low- 
latency quantum storage. 

Work supported by the US ARO under contract No. 
W911NF-11-1-0068, the US NSF under grant No. PHY- 
0903727 (to LV), the IARPA QCS program under con- 
tract No. RC051-S4 (to LV), the IARPA MQCO program 
(to MJB), and the ARC Centre for Engineered Quantum 
Systems, CE110001013. 



ones. In fact, there are \{T s /t) WDD sequences that 
stop at time T s , a very small subset of all 2 Ts / T possible 
digital DD sequences, enabling an otherwise intractable 
numerical BADD search up to T s using Eq. (1). Repre- 
sentative results are shown in Fig. 4, where for each T s 
all WDD sequences with given r are evaluated and those 
with the lowest residual error are selected. As T s grows, 
minimal-error sequences (shown) are found to consist of 
a few concatenation steps (increasing a p of the base se- 
quence) followed by successive repetitions of that fixed 
base sequence. Accordingly, mid-sequence interrupts at 
times corresponding to multiples of T p are permitted in 
low-error long-time storage scenarios without substantial 
performance degradation. In addition to providing a di- 
rect means of finding optimal long-time DD sequences, 
these numerical search results thus support our key ana- 
lytical insights as to use of periodic sequences. 

While our analysis has assumed perfect n rotations, 
various pulse non-idealities are inevitably present in prac- 
tice, their accumulated effect potentially resulting in se- 
vere performance degradation. Beside minimizing the 
number of applied pulses (thus using longer intervals as 
noted), several approaches can be envisioned to miti- 
gate the impact of the dominant pulse imperfections. In 
particular, Euler DD cycles [2] may provide a starting 
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